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Abstract 

Mutualisms between microbes and insects are ubiquitous and facilitate exploitation of various trophic niches by 
host insects. Dictyopterans (mantids, cockroaches and termites) exhibit trophisms that range from omnivory to strict 
wood-feeding and maintain beneficial symbioses with the obligate endosymbiont, Blattabacterium, and/or diverse 
gut microbiomes that include cellulolytic and diazotrophic microbes. While Blattabacterium in omnivorous Periplaneta 
is fully capable of provisioning essential amino acids, in wood-feeding dictyopterans it has lost many genes for their 
biosynthesis {Mastotermes and Cryptocercus) or is completely absent {Heterotermes). The conspicuous functional 
degradation and absence of Blattabacterium in most strict wood-feeding dictyopteran insects suggest that alternative 
means of acquiring nutrients limited in their diet are being employed. A 16S rRNA gene amplicon resequencing 
approach was used to deeply sample the composition and diversity of gut communities in related dictyopteran insects 
to explore the possibility of shifts in symbiont allegiances during termite and cockroach evolution. The gut microbiome 
of Periplaneta, which has a fully functional Blattabacterium, exhibited the greatest within-sample operational taxonomic 
unit (OTU) diversity and abundance variability than those of Mastotermes and Cryptocercus, whose Blattabacterium have 
shrunken genomes and reduced nutrient provisioning capabilities. Heterotermes lacks Blattabacterium and a single 
OTU that was 95% identical to a Bacteroidia-assigned diazotrophic endosymbiont of an anaerobic cellulolytic protist 
termite gut inhabitant samples consistently dominates its gut microbiome. Many host-specific OTUs were identified 
in all host genera, some of which had not been previously detected, indicating that deep sampling by pyrotag 
sequencing has revealed new taxa that remain to be functionally characterized. Further analysis is required to uncover 
how consistently detected taxa in the cockroach and termite gut microbiomes, as well as the total community, 
contribute to host diet choice and impact the fate of Blattabacterium in dictyopterans. 
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Background 

Successful collaborations between microbes and their 
arthropod hosts have enabled exploitation of a wide 
range of trophic niches. Maternally inherited intracellu- 
lar bacteria and/or gut-inhabiting microbial assemblages 
involved in host nutrition are found in insects that per- 
sist on specialized (e.g. plant juices, animal blood, wood 
or decaying plant material) or omnivorous diets (Moran 
et al. 2008; Hongoh 2010). Cockroaches (Dictyoptera: 
Blattaria) are host to both an obligate intracellular sym- 
biont and a species-rich gut microbial community, while 
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most termites (Dictyoptera: Isoptera) host only the latter 
(Lo et al. 2003; Hongoh 2010; Schauer et al. 2012). 
Nearly all cockroaches harbor the intracellular mutualist, 
Blattabacterium sp. (heretofore referred to as Blattabac- 
terium), which has codiversified with their insect hosts 
over millions of years and in most cases they can recycle 
nitrogen from ammonia and urea metabolic wastes into 
the production of essential and nonessential amino acids 
(Sabree et al. 2009). Termites share a common ancestor 
with the wood roach Cryptocercus spp. (heretofore re- 
ferred to as Cryptocercus) via Mastotermes darwiniensis 
(heretofore referred to as Mastotermes), which forms the 
basal branch in termite phylogenies, and as such these 
species represent transitional stages between cockroaches 
and termites. When compared to their modern cockroach 
relatives, termites sport several striking distinctions that 
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include advanced social behaviors, dramatic physiological 
modifications (e.g. enlarged hindguts, perpetual neoteny, 
thin cuticles) (Nalepa 2011) and, with Mastotermes as the 
only known exception, loss of the heritable Blattabacter- 
ium symbiont (Bell et al. 2007). Cryptocercus and Masto- 
termes, also exhibit many of these characteristics and 
retain Blattabacterium (Neef et al. 2011; Sabree et al. 
2012). It is surprising that the genomes of Blattabacter- 
ium in these two host insects lack genes involved in the 
biosynthesis of some essential amino acids, given that 
their hosts are limited for nitrogen given that they thrive 
on wood-based diets typically low in protein (0.03-0.7% 
nitrogen (Merrill & Cowling 1966; Tayasu et al. 1994)). 
Many intracellular bacterial mutualists of phytophagous 
insects with highly reduced genomes still retain essential 
amino acid biosynthesis pathway presumably to supple- 
ment their low nitrogen diet (Sabree et al. 2013). Thus, 
the functional deterioration of Blattabacterium in Crypto- 
cercus and M. darwiniensis and its complete elimination 
in other termites suggest alternative means of obtaining 
and retaining nitrogen in these insects. 

Emergent mutualisms with gut microbes that can pro- 
vide the same nutrients as the ancient endosymbiont 
while conferring new functions may have been a context 
for functional deterioration and eventual loss of Blatta- 
bacterium in Cryptocercus and Mastotermes and ter- 
mites, respectively. Supporting evidence in this regard 
would be the discovery of nitrogen fixation in termites 
(Benemann 1973), its association with hindgut bacteria 
(Yamada et al. 2007; Potrikus & Breznak 1977; Kudo 
et al. 1998; Ohkuma et al. 1999), and identification of 
genes underlying nitrogen fixation and essential amino 
acid biosynthesis in hindgut bacteria in various termite 
species (Wertz et al. 2012; Isanapong et al. 2012). Add- 
itionally, 'lower' termites and Cryptocercus are aided in 
their trophic specialization on wood by lignocellulosic 
hindgut microbes that include both bacteria (Hongoh 
2010; Matteotti et al. 2011; Abt et al. 2012) and protists 
(Tartar et al. 2009; Scharf et al. 2011; Carpenter et al. 
2011; Tamschick & Radek 2013). These protists are 
themselves hosts to intra- and extra-cellular bacterial 
symbionts (Hongoh et al. 2008a; Hongoh et al. 2008b; 
Desai & Brune 2012; Strassert et al. 2012). Since Blatta- 
bacterium can neither degrade cellulose nor fix nitrogen, 
acquisition of organisms capable of these functions en- 
ables their host to exploit an abundant dietary substrate 
and represents a significant improvement on resource 
utilization afforded by Blattabacterium. 

A 16S rRNA gene amplicon resequencing approach 
(Schloss & Handelsman 2003; Tringe & Hugenholtz 
2008) was used to deeply sample the composition and 
diversity of gut communities in dictyopteran insects to 
explore the possibility of shifts in symbiont allegiances 
during the evolution of termites from cockroaches. We 



expected that if Cryptocercus and termites are reliant upon 
their gut microbiomes in the context of a functionally di- 
minished or absent Blattabacterium, respectively, then their 
community profiles will be stable across intraspecific indi- 
viduals and abundant bacterial operational taxonomic unit 
(OTU) therein would be consistently detected. Additionally, 
metabolic profiles and available genomic information from 
taxa related to abundant gut microbiome members de- 
tected in this study were used to make general inferences 
about their possible roles in host trophic ecology. 

Results and discussion 

Sampling effort, community diversity analysis and 
general taxonomic profiling 

The gut microbiota of two termite and two cockroach spe- 
cies were surveyed by pyrosequencing to identify and 
quantify community membership, and to compare intra- 
and inter-host community profiles in the context of 1) the 
host's diet, 2) presence or absence of Blattabacterium and 
3) the inferred functional capacity of Blattabacterium. A 
total of 276,850 high-quality pyrotags representative of 
16S rRNA gene V6-V9 region amplicons were clustered 
into 1,152 OTUs and included in subsequent analyses 
(Additional file LTable SI). Observed and estimated OTU 
richness and rarefaction analyses suggest that near- 
comprehensive sampling of the termite and Cryptocercus 
microbiota could be achieved in fewer than 35,000 pyro- 
tags, but this may not be sufficient for the Periplaneta 
microbiota due to the relatively high OTU richness 
therein (Table 1; Additional file 2: Figure SI). 

Sequences representative of each of the 1,152 OTUs were 
taxonomically classified by comparing them to SILVA (ver- 
sion 108;) and nr (accessed May 6, 2012) databases using 
blastn. OTUs were assigned to over 24 bacterial phyla that 
include Fusobacteria, Deferribacteres, Cyanobacteria, Verru- 
comicrobia, Elusimicrobia, TM7, Planctomycetes, Tenericutes, 
Synergistetes, Actinobacteria, Spirochaetes, Proteobacteria, 
Bacteroidetes, and Firmicutes (Figure 1; Table 2). Single 
OTUs assigned to the Bacteroidia and Clostridia clas- 
ses were abundant (defined as >1%) in every sample, 
and, when combined, predominant in all four host spe- 
cies (Additional file 3: Table S2). Bacilli, Elusimicrobia, 
Deltaproteobacteria and Spirochaetes OTUs were de- 
tected, but they were not uniformly abundant, in every 
sample. With some exceptions, OTUs assigned to many 
of the remaining families were not abundant (i.e. <1%) 
and/or were absent from some samples or sample groups. 

Predominant taxa are host specific and consistently 
abundant 

OTUs assigned to the Bacteroidia predominate in nearly 
all host gut communities but few were shared between 
cockroaches and termites (Table 2). Similarly, Clostridia 
were abundant in every host, but to lesser degree 
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Table 1 Insect gut community alpha diversity 



Sample 


Pyrotags 


Observed OTUs 


Species richness 


OTU obs /OTU ACE 


Species evenness 


Average distance 1 


Ha01 


33427 


189 


191.19 (+/- 4.5) 


0.99 


2.50 (+/- 0.04) 




Ha02 


30441 


150 


171.85 (+/- 17.0) 


0.87 


1.18 (+/- 0.01) 


Ha: 0.23 


Ha03 


34266 


191 


199.90 (+/- 9.6) 


0.96 


1 .85 (+/- 0.02) 




Md01 


16610 


287 


294.43 (+/- 7.7) 


0.97 


1 8.44 (+/- 0.6) 




Md02 


14964 


303 


316.55 (+/- 11.1) 


0.96 


35.33 (+/- 1 .0) 


Md: 0.38 


Md03 


14086 


254 


268.44 (+/- 1 1 .6) 


0.95 


20.41 (+/- 0.5) 




Cp01 


12476 


319 


329.65 (+/- 9.4) 


0.97 


39.65 (+/- 1 .6) 




Cp02 


14482 


321 


338.55 (+/- 1 2.6) 


0.95 


24.35 (+/- 1 .0) 


Cp: 0.31 


Cp03 


12504 


334 


341.44 (+/- 7.5) 


0.98 


34.71 (+/- 1.4) 




PaW01 


20767 


275 


314.39 (+/- 22.7) 


0.87 


1 6.08 (+/- 0.4) 




PaW02 


2081 


225 


282.93 (+/- 29.9) 


0.80 


32.38 (+/- 4.2) 


PaW: 0.57 


PaW03 


30565 


391 


401.13 (+/- 9.2) 


0.97 


26.39 (+/- 0.7) 




PaLOl 


4208 


299 


342.61 (+/- 23.2) 


0.87 


54.68 (+/- 4.5) 




PaL02 


11801 


364 


380.07 (+/- 1 1 .6) 


0.96 


44.60 (+/- 2.0) 


PaL: 0.60 


PaL03 


24172 


384 


400.78 (+/- 1 2.3) 


0.96 


60.96 (+/- 1 .3) 





Species richness was estimated using the abundance coverage estimator (ACE) and evenness was demonstrated as the inverse of Simpson's estimation of 
evenness. 95% confidence indicators are noted in parenthesis for both measures. Ha-Heterotermes aureus, Md-Mastotermes darwiniensis, Cp-Cryptocercus 
punctulatus, Pa\N-Peripianeta americana wild-caught, PaL-P. americana lab-reared. 

1 -Average within-group sample distances were calculated using M u It i- Response Permutation Procedures using the Bray-Curtis distance measure. 
Test statistic: -8.719, observed delta: 0.4142, expected delta: 0.8829, chance-corrected within group agreement: 0.5309, P-value: 7.0e-8. 

(representing <4% of the total OTUs) in H, aureus. Like yet few OTUs were shared, the observed host-defined 

the Bacteroidia-assigned OTUs, none were shared across sample clustering in the NMS analysis is not unexpected 

both termites and cockroaches. Given that 59-96% of the (Figure 2). It is not surprising to find members of both 

pyrotags in each sample could be assigned to either class, classes in all of the gut microbiomes of hosts in this study 
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Figure 1 Bacteroidia and Clostridia predominate in most cockroach and termite gut microbiota. Heatmap depicts the relative abundances 
of pyrotags assigned to each taxonomic class; total number of pyrotags representative of each sample are parenthesized. 



Table 2 Relative abundance of bacterial families present within insect gut communities 



Taxonomic Ha01 Ha02 Ha03 Md01 Md02 Md03 Cp01 Cp02 Cp03 PaL01 PaL02 PaL03 PaW01 PaW02 PaW03 
designation 

Actinobacteria; 1.06(8) 0.23(7) 1.01(8) 2.74(8) 2.74(10) 4.14(9) 0.70(10) 0.37(10) 0.49(9) n.d. (0) n.d. (0) n.d. (1) 0.18(6) 0.91(7) 0.51(8) 
Actinobacteria 

Actinobacteria; 0.05 (2) 0.02(2) 0.01(2) 0.55 (3) 1.21 (3) 0.47(3) 1.00 (4) 1.07(4) 0.86 (4) 0.12 (2) 0.08 (2) 0.13(2) 4.49 (2) 1.44 (2) 0.99 (2) 
Coriobacteria 

Bacteroidetes; 75.83 (21) 95.68(17) 83.53 (18) 42.41(52) 28.31 (51) 33.95(50) 36.66 (46) 48.57(50) 44.71(50) 41.02 (89) 58.21(88) 56.59(100) 1.51 (50) 19.75 (55) 14.83 (88) 
Bacteroidia 

Candidate 0.04(2) 0.01(1) 0.04 (2) 1.53 (13) 2.22(14) 0.94(11) 1.08 (12) 1.09(10) 0.65(8) 0.26 (4) 0.09 (3) 0.31(6) 0.31 (5) 1.87 (7) 0.29 (6) 
division TM7 

Deferribacteres; n.d. (0) n.d. (0) n.d. (0) n.d. (0) n.d. (0) n.d. (0) 3.51 (1) 4.81(1) 1.45(1) 0.02 (1) 0.23 (2) 0.23(2) n.d. (0) n.d. (0) 0.05(2) 
Deferribacteres 

Elusimicrobia; 0.04(2) 0.01(2) 0.02 (2) 4.31 (8) 1.87 (9) 6.07(7) 0.57 (2) 0.26(3) 0.8 (2) 2.16 (5) 1.07 (5) 0.26(3) 0.01 (1) 0.05 (1) 0.1(3) 
Elusimicrobia 

Firmicutes; 5.92(9) 0.2(7) 0.81(11) 4.18 (5) 3.47 (5) 1.54 (4) 5.79 (9) 7.88(8) 7.98(8) 20.29 (7) 0.24 (5) 0.01 (2) 1.54 (8) 4.13 (6) 24.8 (10) 
Bacilli 

Firmicutes; 4.00(67) 1.05(50) 3.11(67) 24.75 (97) 38.2 (109) 17.7 (84) 35.02(140) 24.5 (138) 29.19 (147) 22.98 (117) 22.4 (151) 15.86(151) 61.81(142) 48.82 (100) 44.06 (186) 
Clostridia 

Firmicutes; n.d. (0) n.d. (0) n.d. (0) 0.16 (2) 0.43 (2) 0.07 (2) 0.05 (1) 0.09(1) 0.02 (1) 0.9 (6) 0.23 (7) 0.55(7) 1.58 (5) 0.96 (3) 0.4 (5) 
Erysipelotrichi 

Fusobacteria; n.d. (0) n.d. (0) n.d. (0) 3.00 (1) 4.25 (1) 6.12 (1) n.d. (0) n.d. (0) n.d. (0) 0.02(1) n.d. (0) n.d. (0) 13.71(1) 1.11(1) 0.22(1) 
Fusobacteria 

Lentisphaerae; n.d. (0) n.d. (0) n.d. (0) 0.30 (3) 0.25 (3) 0.1 (2) n.d. (0) n.d. (0) n.d. (0) 0.26 (4) 1.49 (5) 3.17 (5) 0.18 (2) n.d. (0) 0.08(2) 
Lentisphaeria 

Planctomycetes; 0.43 (2) 0.01(1) 0.13 (2) 0.24 (3) 0.12 (2) 0.19 (3) 0.24 (1) 0.26(2) 0.12 (2) 0.12 (3) 0.29 (6) 0.14(5) 0.83 (5) 0.19 (4) 1.85 (6) 
Planctomycetacia 

Planctomycetes; 0.04(1) 0.03(1) 0.11(1) n.d. (0) n.d. (0) n.d. (0) 2.65 (4) 1.49(3) 2.36(4) 0.05 (2) 0.28 (5) 0.67(5) 3.29 (4) 1.06 (4) 2.65 (5) 
vadinHA49 

Proteobacteria; 0.08(4) 0.06(3) 0.04(4) 0.14(6) 0.61(9) 0.24(4) 1.62(15) 1.17(14) 0.86(14) 1.88(9) 2.3(14) 1.42(16) 0.03(3) 0.1(1) 0.04(5) 

Alphaproteo- 

bacteria 

Proteobacteria; 0.70(5) 0.42(6) 0.71(6) 0.73 (5) 1.22 (5) 0.6(4) 1.52 (7) 1.36(9) 2.39(9) 0.71 (3) 0.36 (6) 0.19(6) 3.39 (6) 0.43 (1) 0.27 (3) 
Betaproteobacteria 

Proteobacteria; 0.35(5) 0.06(5) 0.41(6) 0.02 (1) 0.03 (1) 0.15 (1) 0.85 (9) 0.16(7) 0.46 (8) 3.26 (20) 5.94 (26) 10.38 (32) 5.84 (20) 12.88 (15) 5.69 (26) 

Deltaproteo- 

bacteria 

Proteobacteria; 0.59(2) 0.01(1) 0.1(2) 0.22 (2) 0.34 (3) 0.23 (2) 0.35 (4) 0.19(6) 0.42 (6) 0.64 (2) 0.18 (1) n.d. (0) 0.03(2) 0.05(1) 0.04(2) 

Gammaproteo- 

bacteria 

Spirochaetes; 8.84 (39) 1.87 (32) 8.34(40) 9.67 (45) 5.93(39) 20.02(42) 4.74(15) 4.47(15) 4.88(16) 2.38 (7) 2.69 (11) 2.96(12) 0.43 (1) 1.3 (4) 0.72 (6) 
Spirochaetes 



Table 2 Relative abundance of bacterial families present within insect gut communities (Continued) 



Synergistetes; 
Synergistia 


1.24 (5) 


0.12 (5) 


0.75 (5) 


0.55 (7) 


1 .09 (8) 


0.6 (7) 


1.31 (5) 


0.77 (5) 


0.79 (5) 


0.26 (2) 


0.47 (4) 


1.21 (4) 


0.35 (3) 


1 .78 (2) 


0.81 (3) 


Teneticutes; 
Mollicutes 


0.01 (1) 


0.02 (1) 


0.02 (1) 


1 42 (7) 


5.02 (7) 


4.87 (5) 


0.66 (7) 


0.33 (6) 


0.47 (8) 


0.64 (4) 


1 .55 (4) 


1 -25 (5) 


0.22 (3) 


0.86 (3) 


0.25 (5) 


Verrucomicrobia; 
Opitutae 


0.02 (1) 


0.01 (1) 


0.03 (1) 


2.30 (5) 


1.3 (4) 


1.27 (5) 


0.55 (6) 


0.53 (7) 


0.31 (8) 


0.07 (1) 


0.08 (2) 


0.05 (2) 


n.d. (0) 


0.05 (1) 


0.04 (1) 


Other bacterial 
families 
(21 families) 


0.77 (8) 


0.19 (5) 


0.83 (8) 


0.52 (10) 


1.14 (14) 


0.19 (5) 


0.98 (12) 


0.38 (12) 


0.61 (14) 


1 .69 (9) 


1.2 (14) 


3.88 (16) 


0.2 (4) 


2.07 (5) 


1.22 (12) 


Archaea 
(4 families) 


n.d. (0) 


n.d. (0) 


n.d. (0) 


0.25 (3) 


0.25 (3) 


0.54 (2) 


0.14 (2) 


0.26 (3) 


0.17 (3) 


0.26 (1) 


0.64 (3) 


0.73 (2) 


0.06 (2) 


0.19 (2) 


0.1 (4) 



Values represent the percentage of all taxonomically-assigned OTUs within each sample. Parenthesized values represent the real number of OTUs within each sample assigned to that taxonomic group, 
n.d.-not detected. 
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termites cockroaches 

o Heterotermes aureus Periplaneta americana, lab-reared 

Mastotermes darwinensis v Periplaneta americana, wild-caught 
A Cryptocercus punctulatus 

Figure 2 Insect community samples form host-defined clusters. Nonmetric multidimensional scaling was used to visualize community 
diversity and relatedness. A randomly subsampled, size-normalized (n = 2,050) pyrotag abundance data table was used. 



as cultivated bacteroidia and Clostridia 1) range from aero- 
tolerant to strict anaerobes, exhibit a wide range of poten- 
tially host-beneficial carbon fermentative metabolisms, 
including cellulose degradation by various clostridial gen- 
era (Tracy et al. 2012), 2) are well-adapted to the dy- 
namic gut environment and 3) are common amongst 
the gut flora of healthy invertebrates and vertebrates 
(Schleifer 2009; Engel & Moran 2013). 

OTUs assigned to the Tenericutes and Elusimicrobia 
were abundant only in Mastotermes. Members of the 
Tenericutes, namely Phytoplasma, Mycoplasma and Spir- 
oplasma are well-adapted to reside in various animal 
and plant hosts largely as pathogens (Gamier et al. 2001; 
Gasparich 2010) but a recent report of Spiroplasma con- 
ferring protection against parasitic nematodes in some 
Drosophila species (Haselkorn et al. 2013) indicates other 
interactions between tenericutes and their eukaryotic 
hosts are possible. 

Members of the Elusimicrobia were initially detected 
in Reticulitermes speratus hindguts (Ohkuma & Kudo 
1996; Hongoh et al. 2003) and have since been detected 
in the guts of various phytophagous insects and 



terrestrial habitats (Herlemann et al. 2007). Elusimicro- 
bium minutum was the first cultivated member of this 
clade (Geissinger et al. 2009) and this strict anaerobe 
generates ATP through typical fermentative pathways 
and generates some amino acids but likely relies upon 
other microbiome members, the host or the host's diet 
for additional required nutrients. In contrast, production 
of many essential amino acids using ammonia by an Elu- 
simicrobia endosymbiont of a cellulolytic protist inhabit- 
ant of R. speratus was inferred from an analysis of its 
genome and they are accessible to the host via digestion 
of the protist or by uncharacterized transport mecha- 
nisms (Hongoh et al. 2008a). 

Although Spirochaetes have been well documented in ter- 
mite gut communities (Paster & Dewhirst 2000; Rosenthal 
et al. 2011; Kohler et al. 2012), no single OTU was abun- 
dant in either termite host species examined. Hydrogen 
metabolism has been proposed as one of the functions of 
spirochaetes in termite guts (Rosenthal et al. 2011) but fur- 
ther characterization is necessary. Combined OTUs 
assigned to the Spirochaetes were present and >1% of 
the sequence reads in all insects except two wild-caught 
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P. americana samples, yet only three occurrences were 
noted where single OTUs were >1% of the sequence reads 
for a sample (Additional file 3: Table S2). Since the OTUs 
were defined at >95% sequence identity, the Spirochaetes 
detected in the termites exhibited greater within-sample 
diversity and evenness of abundance than that observed 
for OTUs assigned to the Bacteroidia or Clostridia classes 
in the same termites. 

Protist-associated OTUs predominate the microbiota of 
insects lacking fully functional Blattabacterium 
endosymbionts 

Nearly all cockroaches harbor the obligate intracellular 
mutualist, Blattabacterium, and its absence in all ter- 
mites except Mastotermes presents an opportunity to ex- 
plore the possibility of mutualist replacement within the 
dictyopterans. We hypothesized that loss of Blattabac- 
terium could only be tolerated if newly acquired mi- 
crobes 1) were abundant in the gut (a proxy for being 
well-adapted to the host environment and an important 
member of the community), 2) could be reliably trans- 
mitted to offspring and 3) capable of supplying nutri- 
ents required by the host. Results from this survey of 
the microbiota indicate that microbiota samples from 
Heterotermes, which lacks Blattabacterium, had the 
least amount of OTU diversity and evenness and exhib- 
ited the least community profile variability, as indicated 
by the low average distance (Table 1). Conversely, 
microbiota samples from P. americana, which harbors 
Blattabacterium with all of its essential amino acid bio- 
synthesis pathways intact, were more variable for rela- 
tive abundance, OTU richness, evenness and thus 
exhibited the greatest distance between samples from 
the same host. Average distances for the microbiota 
samples from Cryptocercus and Mastotermes, whose 
Blattabacterium associates are deprived of many of the 
aforementioned biosynthetic pathways, were intermedi- 
ate to those of Heterotermes and Periplaneta. 

Accounting for a large proportion of the Heterotermes 
microbiota is a single Bacteroidia-assigned OTU (OTU_ 
0035) that had a best hit (95% identity) to a nitrogen- 
fixing (diazotrophic) endosymbiont of the cellulolytic 
trichonymphid protist, Pseudotrichonympha grassii, that 
resides in Coptotermes formosanus hindguts (Hongoh 
et al. 2008b) (Additional file 3: Table S2). The prevalence 
of OTU_0035 in the H. aureus samples suggests that it 
may be an essential member of the community, which 
would be reasonable if it too is diazotrophic and partnered 
with a cellulolytic protist. Soldiers and young instars ob- 
tain nutrients and gut microbes from the hindgut fluids of 
conspecifics via proctodeal trophallaxis (Machida et al. 
2001; Nalepa et al. 2001). This practice would ensure 
the reliable communication of the gut microbiota that in- 
clude oxygen-sensitive microbes, some of whom can fix 



nitrogen, a process that is inactivated by oxygen. The tight 
mutualism between the N 2 -fixing endosymbiont and the 
cellulolytic P. grasii combines functions essential for the 
termite to thrive on plant-based carbon sources in a single 
trophic mutualism. If the role of Blattabacterium is to 
provision amino acids and vitamins missing in the host's 
diet, its role would be tangential in the presence of the cel- 
lulolytic protist-diazotrophic bacteria symbiosis because 
the N2-fixing protist endosymbiont can generate the same 
assortment of nutrients and Blattabacterium cannot de- 
grade cellulose or fix atmospheric nitrogen. 

Both Mastotermes and Cryptocercus have cellulose- 
based diets, exhibit trophallaxic behavior and still retain 
Blattabacterium, albeit with reduced genomes and di- 
minished nutrient provisioning abilities. Unlike Blatta- 
bacterium in other cockroaches, which are capable of 
making all ten essential amino acids, the endosymbiont 
in Mastotermes and Cryptocercus encode the genes for 
the production of only five or six essential amino acids, 
respectively (Neef et al. 2011; Sabree et al. 2012). In most 
cases members of the Bacteroidia were the most abundant 
in Cryptocercus and Mastotermes, but no single OTU was 
as abundant as OTU_0035 in Heterotermes. 

Two abundant, protist-associated bacteroidial OTUs 
(OTU_1551 and OTU_0214) in Mastotermes were only 
93% identical to OTU_0035, which suggests little gene 
flow between these bacteroidial taxa and isolation of 
their protist hosts within their respective termite hosts 
would support this. The devescovinid Mixotricha para- 
doxa is abundant in Mastotermes but the functions of its 
associated bacteroidia remain to be determined. Add- 
itional OTUs detected in Mastotermes that were related 
to bacterial symbionts of gut protists were taxonomically 
assigned to the Elusimicrobia and Tenericutes. Protists 
are abundant symbiotic inhabitants of termites (Ohkuma 
& Brune 2011; Desai & Brune 2012), which likely con- 
tributes to the relative abundance of their associated 
bacteria. The low percent identity to available taxon- 
assigned sequences suggests that this is the first pub- 
lished detection of the TeMen'cKtes-assigned gut protist 
symbiont in M. darwiniensis. It is possible that the M. 
darwiniensis Elusimicrobia detected in this study plays a 
role in life of its host that is analogous to that of R. sper- 
atus protist symbionts given that their hosts have similar 
diets. 

Many of the Bacteroidia and Clostridia-assigned OTUs 
in the P. americana gut communities had best hits to 
amplicons obtained from either Shelfordella lateralis 
(Schauer et al. 2012) or various termites. S. lateralis and P. 
americana are both part of the Blattidae (Blattoidea) fam- 
ily, which suggests that the shared OTUs may represent 
taxa that are well-adapted to their cockroach hosts and 
were acquired following divergence of the Shelfordella- 
Periplaneta (Blattoidea) and Cryptocercus (Blaberoidea) 
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clades. Deep sampling of gut microfloral diversity of other 
host taxa in both of these superfamilies is necessary to test 
this hypothesis. The role of taxa shared by Periplaneta and 
Shelfordella in cockroach host development and/or troph- 
isms remain to be characterized. Finally, the relative abun- 
dance of pyrotags representing the abundant OTUs 
assigned to the Bacteroidia and Clostridia classes varied 
significantly (T-test: p < 0.0003) between the samples ob- 
tained from wild-caught and lab-reared P. americana. This 
supports the hypothesis that P. americana harbors some 
host-specific bacteria, as evidenced by the ordination ana- 
lyses, but diet and/or habitat can impact their prevalence. 

Ammonia-oxidizing and sulfate-reducing bacteria 
in cockroaches 

OTUs assigned to the Planctomycetes and Deltaproteobac- 
teria were largely shared by and unique to the cockroaches. 
Known physiological characteristics of Planctomycetes 
members are the lack of peptidoglycan in their cell walls, 
intracellular compartmentalization, sterol biosynthesis, and 
budding reproduction (for review (Fuerst & Sagulenko 
2011)). Additionally, so-called 'annamox' members are 
chemoautotrophic anaerobes capable of oxidizing ammo- 
nia to dinitrogen and are being successfully exploited for 
energy-efficient removal of nitrogenous wastes in wastewa- 
ter treatment (Siegrist et al. 2008; Shi et al. 2013). Ammo- 
nia comprises much of the nitrogenous wastes externally 
excreted by cockroaches (Cochran 1985) and annamox 
Planctomycetes present in their hindguts may utilize this 
surplus to generate ATP via the annamoxosome (van 
Niftrik & Jetten 2012) and, by effect, may have a detoxify- 
ing effect by reducing the concentration of ammonia 
present in the hindgut. 

Deltaproteobacteria-assigned OTUs in Cryptocercus 
and Periplaneta hindguts had best hits to amplicons that 
were assigned to the Desulfobacteriaceae and Desulfovi- 
brionaceae families. Many cultivated members of both 
groups are strict anaerobes residing in marine sedi- 
ments that are capable of coupling sulfate-reduction 
and energy production (Muyzer & Stams 2008). Addition- 
ally, Candidatus "Desulfovibrio trychonymphae", an endo- 
symbiont of an anaerobic R. speratus-mhabiting protist, 
has also been shown to have and express genes involved 
in sulfate-reduction (Sato et al. 2009), but it is not clear if 
this is their primary function in the hindgut microbiome. 

Conclusions 

The gut microbiomes of dictyopteran insects surveyed in 
this study are comprised of many bacteria from the same 
class-level taxonomic groups (e.g. Bacteroidia, Clostridia, 
Spirochaetes and Bacilli) but distinct sub-lineages were ob- 
served when OTUs were resolved at the >95% sequence 
identity cutoff, indicating the presence of many relatively 
abundant host-specific taxa. The gut communities of the 



strict wood-feeding insects had bacterial taxa known to be 
associated with cellulolytic protists and fewer shared OTUs 
while the OTU diversity was generally greater and more 
variable in the omnivorous P. americana. If the abundant 
endosymbiont of the H. aureus cellulolytic trichonymphid 
protist is diazotrophic and capable of provisioning amino 
acids to its protist host that is itself digested by the termite, 
then acquisition of this and other gut bacterial-protist or 
bacterial symbioses that facilitated wood-feeding could 
have contributed to the shift in nutritional reliance from 
Blattabacterium to the gut microbiome. M. darwiniensis 
and Cryptocercus have microbiomes that are more diverse 
than that of H. aureus, harbor cellulolytic protists with bac- 
terial symbionts and also sport many of the host physio- 
logical and behavioral modifications observed in H. aureus. 
Given that M. darwiniensis and Cryptocercus are sister 
taxa, both harbor Blattabacterium spp. that are function- 
ally reduced and M. darwiniensis is basal to termites, it is 
possible that these host gut microbiomes represent inter- 
mediate stages of a more stable community that is essential 
for wood-feeding. Unlike in M. darwiniensis and Crypto- 
cercus, the P. americana Blattabacterium is equipped to 
provision vitamins and a near- complete suite of amino 
acids to its host, which may reduce its reliance upon the 
gut microbiome for these functions. Elevated average dis- 
tances between within-group samples for P. americana gut 
microbiota membership suggest this possibility but the 
presence of shared OTUs between the wild-caught and 
lab-reared P. americana indicates a host-specific micro- 
biota that is present regardless of diet or lifestyle. The func- 
tions of these members remain to be determined. 

Deep sequencing has helped to identify a number of 
previously undetected taxa in this study, and others 
seeking to profile cockroach and termite microbiomes 
(Kohler et al. 2012; Schauer et al. 2012; Huang et al. 2013; 
Boucias et al. 2013) that likely represent novel strains, 
species or genera, indicating that the guts of dictyopteran 
insect, of which there are about 8,500 species, contain a 
wealth of novel bacterial diversity. It is clear that further 
functional characterization of abundant cockroach and 
termite gut microbiome members is necessary and will 
likely reveal some new biological activities. 

Methods 

DNA extraction and multiplexed sample preparation 

Entire guts were dissected from fresh or ethanol-preserved 
specimens from the following sources: Heterotermes aureus 
(Tucson, Arizona, USA; July 2008), Mastotermes darwi- 
niensis (Marlow Lagoon, Northern Territory, Australia), 
Cryptocercus punctulatus (Mountain Lake, Virginia, USA), 
wild-caught Periplaneta americana (near the University of 
Arizona gymnasium, Tucson, Arizona, USA, July 2010, 
33°C and 25% relative humidity), lab-reared P. americana 
(lab colony fed on dog food containing 28% amino acids, 
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provided water ad libatum and maintained at 24C and 
35% relative humidity). Three adult individuals captured 
from each habitat were sampled. DNA was prepared using 
the Power Soil DNA Isolation Kit (MoBio, San Diego, 
USA) according to the supplied protocol. 50 ng of template 
DNA was used in PCR amplifications performed in 
triplicate in 30 uL reactions containing 0.4 uM bacteria- 
specific forward primer (TAXX-926 F: 5'-{adapter}- 
{barcode}-AAACTYAAAKGAATTGACGG-3 '; (Lane 1991; 
Engelbrektson et al. 2010)) that were uniquely barcoded 
for multiplexing samples, 0.4 uM nonbarcoded universal 
primer (TB-1392R: 5 '-{adapterJ-TACGGYTACCTTGT 
TACGACTT-3'; (Ferris et al. 1996; Engelbrektson et al. 
2010)) (see Additional file 4: Materials for primer se- 
quences), 0.55 U Phusion Taq DNA polymerase (New 
England Biolabs, Massachusetts, USA), and 1 mM dNTP 
mix (Promega, Wisconsin, USA). Selected primers were 
used to amplify the hypervariable V6-V9 regions of the 
16S rRNA gene (Sogin et al. 2006; Roesch et al. 2007). 
Reactions were initially denatured at 98°C for 1 min, 
followed by 25 cycles of 98°C for 10 sec, 55°C for 
10 sec and 72°C for 15 sec and a final, single cycle of 
72°C for 10 min. Amplification was confirmed by gel 
electrophoresis, and PCR products from amplifications 
performed in triplicate were pooled, purified using 
AmPure magnetic beads (Beckman Coulter, Indianapolis, 
USA) and quantified using the Qubit fiuorometer (Life 
Technologies, New York, USA). Barcoded PCR products 
were combined, at a final concentration of 0.45 ng per re- 
action, into a single, multiplexed sample that was submit- 
ted to the University of Arizona Genomics Center for 
pyrosequencing on a Roche 454 FLX-Titanium system. 

Pyrotag processing and analysis 

402,054 raw, barcoded amplicon sequences ("pyrotags") 
were obtained from the pyrosequencing run. These were 
processed within the CLC Genomics Workbench (www. 
clcbio.com) to trim low-quality regions (<27 Phred score) 
from each pyrotag and to remove reads that were less 
than 400 bp in length and/or had one or more errors in 
the forward primer or barcode regions. MOTHUR was 
used for further pyrotag processing (version 1.24, 
(Schloss et al. 2009)) (for details, see Additional file 4: 
Materials), and pyrotags having >95% identity were 
clustered into operational taxonomic units (OTUs) in 
MOTHUR (cluster.split, method = furthest). OTUs that 
did not have at least two pyrotags present in at least 
two samples were excluded from our dataset to minimize 
the impact of possible contaminants and extremely rare 
OTUs on our analyses. OTUs were taxonomically classified 
by blastn-based (BLAST+, version 2.2.26; parameters -task 
blastn -outfmt 6 -evalue le-50; (Camacho et al. 2009)) 
searches of Silva (version 108, (Pruesse et al. 2007)) 
and NCBI 'nt' nucleotide databases with representative 



sequences of each OTU. Acceptable alignments included 
only those for which >85% of the query and 'hit' (or 
'subject' in BLAST+ parlance) sequences were aligned. 
Available taxonomic information for the top hits were 
used to define OTUs. OTUs having hits to plastids or 
mitochondria were removed from the analysis. We used 
the abundance-based coverage estimator (ACE) (Chao 
et al. 1993) to conservatively predict the number of 
OTUs in each sample (richness) and the inverse 
Simpson diversity index (Simpson 1949) to estimate 
OTU evenness. Prior to our comparative community 
analyses, we generated a size-standardized subset of the 
original data that, in terms of the relative abundances 
of OTUs in each sample, was not significantly different 
(paired t-test, p > 0.05) from the original dataset by ran- 
domly sampling pyrotags from each OTU to 2,050 per 
sample for all samples. The subsampled dataset was 
used as input to test and visualize within-group and 
between-group differences using Multi-Response Per- 
mutation Procedures (MRPP) (Mielke 1984) and Non- 
metric Multidimensional Scaling (NMS) (Kruskal 1964) 
within the PC-ORD (version 4.25) software package 
(McCune & Mefford 1999). The Bray-Curtis distance 
measure was used in both comparative analyses and a 
combination of low stress and maximum stability was 
sought for the NMS solution. 

Sequencing data accession number 

Representative OTU names and corresponding NCBI 
GenBank accession numbers can be found in Additional 
file 4. 

Additional files 

r *i 

Additional file 1: Table SI. Pyrotag processing data. Numbers indicate 
quantities of reads before and after trimming low quality base calls and 
removal of undersized reads, a-percentage of total remaining reads, 
b-percentage of total remaining OTUs. 

Additional file 2: Figure SI. Insect community sampling analysis. 
Rarefaction curves reflect sampling-without-replacement. A: Ha-Heterotermes 
aureus, B: Md-/Wa5fofermes darwiniensis, C: Pa\N-Periplaneta americana 
wild-caught, D: PaL-P. americana lab-reared, E: Cp-Cryptocercus punctulatus. 
CI- 95% confidence intervals. Number of pyrotags are indicated on the 
x-axis and number of OTUs are indicated on the y-axis. 

Additional file 3: Table S2. Best GenBank matches to abundant OTUs. 
OTUs that were present in all three samples for at least one host sample 
group and represented >1% of the total pyrotags in at least two samples 
per sample group are described. Asterisks indicate that fewer than 1% of 
the total pyrotags for that sample clustered with the corresponding OTU. 
Sequences representative of selected OTUs were compared to a local 'nt' 
database using blastn (parameters: -task blastn -outfmt 6) and sequences 
with the greatest 'max identity' and 'query coverage' were reported. All 
hits were to uncultured bacteria detected by PCR unless otherwise noted. 
Description information was obtained from the GenBank entry 'definition' 
and 'accession' entries. Accession numbers for best hits are parenthesized. 
Habitat information was obtained from the GenBank entry 'source' 
section.%ID- percent identity. A - According to searches of 'nt' and 'silva' 
databases, GenBank entry AM422253.1 is incorrectly annotated as an 
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"Uncultured Eubacterium sp". and would not be a member of the 
Firmicutes. &- Best blast hit was to a bacterial isolate. 

Additional file 4. Representative OTUs. Text document includes 
representative OTU names and corresponding NCBI GenBank accession 
numbers. 
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